Thermo-Plasma Polariton within Scaling Theory of Single-Layer Graphene 
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Electrodynamics of single-layer graphene is studied in the scaling regime. At any finite temper- 
ature, there is a weakly damped collective thermo-plasma polariton mode whose dispersion and 
wavelength dependent damping is determined analytically. The electric and magnetic fields associ- 
ated with this mode decay exponentially in the direction perpendicular to the graphene layer, but 
unlike the surface plasma polariton modes of metals, the decay length and the mode frequency are 
strongly temperature dependent. This may lead to new ways of generation and manipulation of 
these modes. 



There has been renewed interest in the electronic prop- 
erties of the single layer graphene, spurred largely by 
the recent progress on the experimental front yl |2|. 
Most of the experimental [H Q. Q and theoretical study 
[HI3 has focused on the electronic properties in 

the quantum Hall regime, where the massless relativistic 
dispersion of the electrons at low energies leads to effects 
absent in the conventional two dimensional electron gas. 

On the other hand, the electronic properties of the sin- 
gle layer graphene in the absence of the external magnetic 
field are arguably as interesting as in the quantum Hall 
regime. The undoped system has effectively one particle 
per site, but two per unit cell of a honeycomb lattice, 
and in the non-interacting limit the valence and conduc- 
tion bands touch at two inequivalent K points (see e.g. 
[ToL Hlj p. To a very good approximation, the dispersion 
near these points is £k = i>j?|k| up to e ~ leV, and the low 
energy excitations behave as 2D massless Dirac particles. 

In the weak coupling regime, the role of the electron- 
electron (e — e) interactions at zero temperature has been 
studied previously in a scries of papers 0D3G3, but 
because the results are scattered across the literature, it 
is useful to review and interpret the main conclusions 
here. Since there are no mobile carriers at T = 0, the 
Coulomb interactions are not screened and the usual 
Fermi liquid arguments do not guarantee sharp quasi- 
particle excitations at the Fermi point. The strength of 
the e — e interactions can be parameterized by a dimc- 
sionless coupling constant ap — e 2 /(eohvp), where eo is 
the bare dielectric constant due to the polarizability of 
the core electron bands, cyf is of order one. Standard 
renormalization group (RG) arguments show that in 2D 
the charge e is exactly marginal [lil fl6| , and therefore 
any RG flow of ap- is due to the renormalization of Vp. 
In weak coupling the scale dependent vf grows under RG 
as Vp(k) = Vp^ + |^ log [4^1 and therefore ap decreases, 
resulting in the effective weakening of the e — e interac- 
tions. The growth of vp is eventually bounded by c, the 
speed of light, and at the eventual infra-red fixed point 
the full retarded interaction between electrons mediated 
by photons must be taken into account. This fixed point 
is charged, with electrons moving at the speed of light, 



and it is almost certainly not a Fermi liquid. Unfortu- 
nately, its character, while fascinating, remains only of 
academic interest since c/vp ~ 300 and it would be im- 
possible to probe the system across scales differing by 
hundreds of orders of magnitude. The physical proper- 
ties of the system are then given by a large crossover 
regime dominated by the non-retarded Coulomb interac- 
tions, which, at least in weak coupling can be described 
within perturbative RG. 

The absence of the screening at T = leads to the ab- 
sence of any well defined plasma oscillations. However, as 
shown below, at any finite temperature, there is a finite 
non-activated density of thermally excited quasiparticles 
which can screen the Coulomb interactions. Their cou- 
pling to the three dimensional electromagnetic radiation 
leads to a thermo-plasma polariton mode with frequency 
vanishing as w p = c|k|, but crossing over into lo p — 
(k B T/h)Vt{t+g)/y/i + 2g, where g = e 2 iVln2/(27r^ F ), 
N = 4 is the number of two component Dirac species, 
and t — Hvp\k.\/(kBT). The crossover wavelength is in- 
versely proportional to temperature and is of the order 
of a millimeter at room temperature. The electromagen- 
tic fields are confined to the two-dimensional graphene 
sheet within the attenuation length £ a , which prevents 
the energy from radiating outwards. The analytical ex- 
pressions for u>p, £ a , and the propagation length, 5£, are 
the main results of this paper. While surface plasma 
polaritons have been discussed in metals long time ago 
|17| . the strong temperature dependence of the disper- 
sion and the propagation length of the thermo-plasma 
polaritons, as well as their confinement to the atomically 
thin graphene sheet, may open new opportunities in the 
emerging field of plasmonics 0, . 

To justify the above claims, I start by solving the 
Maxwell equations in the presence of a graphene sheet. 
The speed of light c is taken as a variable and not nec- 
essarily equal to the speed of light in the vacuum. This 
allows for situations where the graphene sheet is posi- 
tioned between two dielectrics. In the scaling regime, all 
energy scales can be measured in units of k B T and all 
lengths can be measured in the units of thermal length 
£t = hvp / (ksT), where Vp is the Fermi velocity at the 
nodal points. As mentioned above, the strength of e-e 



interactions is given by the dimensionless coupling con- 
stant, ctF which is the graphenc analog of the fine struc- 
ture constant. Since the logarithmic corrections to Vf do 
not change the main physics discussed below, for the pur- 
poses of this article these corrections will be neglected. 

Within the linear response theory the effective La- 
grangian for the electromagnetic gauge field is 
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where, A^\ z= q — J ^-A tl (k z ) and for the sake of brevity, 
the dependence of A^ on u> and k was suppressed. In the 
Lorenz gau ge the Fourier components of the polarization 
tensors l2fj| are 
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where [i runs from to 3 while a, (3 — 0,1,2; q a — 



{— iu, t»Fk), so q = —uj 



P k . The diagonal ten- 



sor -500 = gn = 522 = 533 = 1 while ig£ = ^g n = 

£s£ = l,sS = 0. 

In the scaling regime, the polarization functions de- 
pend only on u>, vpk, and T. All the energy (length) 
scales can be measured in ksT (£t) since the response 
is due to the excitations around the nodal points, and 
therefore the polarization has the scaling form 



n Q ,6 = U atb (ui,v F k,T) = ij^Va,, 



\k\£ T ;a F ,N ,(5) 



where the dependence on the dimensionless coupling con- 
stant cxf — e 2 / (eohv p) , and the number of nodal points 



N, was included explicitly on the right hand side. So 
far, the only approximation used is the linearity of the 
electromagnetic response and the existence of the scaling 
regime with vanishing anomalous dimension. 

The equations of motion (Maxwell's equations) for the 
gauge field 
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are invariant under rotations inside the plane as well as si- 
multaneous scale transformation of space-time and ther- 
mal length. Therefore in the scaling regime, any solution 
must satisfy 



w(k, k z ,T) 



w(e T \k\,e T k z ). 



(6) 



Clearly the z-component of A^ decouples from the rest 
and obeys the standard wave-equation. To proceed, we 
find the common (normalized) eigenvectors of and 
b pu and decouple the Eq. JfJJl. This orthogonal triad 

is eH, = ei a) = (S 0a - ujn a / k 2 )\k\/ {c\\l\), and 

e { a = eoa/33W|k|, where k m = (w,ck), and e M „ A p is the 
completely antisymmetric rank 4 tensor. Projecting A^ 
along the four eigenvectors we have 
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The Maxwell's equations then decouple to 
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If c 2 k 2 7^ lu 2 and if we let A\\ — J 4 Q 2wc|k|/(c 2 k 2 — uj 2 ), Lorenz condition ~^Aq + kiAi + k z A z = is satisfied 
then the first equation is automatically satisfied and the provided that A z (z) — constant. Fourier transforming 



into the real space along z-direction we finally have 
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The above equations resemble a Schrodinger equation 
for a particle moving along a line and experiencing a 
5— function potential. If the potential is attractive, we 
know that there is one bound state, otherwise there are 
only scattering states. In this case, the potential depends 
on the frequency u, on the inplane wavevector k and on 
the temperature. The bound state solutions correspond 
to the evanescent electromagnetic waves whose vector po- 
tentials A a and A b have the form 

A a . b (Lu k (T), k, z; ^) = A a , b (0) e -M/( 2 ^) (15) 
c 

and whose inverse attenuation lengths satisfy 

II n and 
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The condition -Re [S, a ,b] > thus specifies a region in 
(w,k)-space where the curve given by the solution of 
Eq.jnj must lie. 

Using the scaling form of the polarization functions, 
Eq.©, gives the scaling form of the attenuation lengths 

Ub(u,v F k,T)=t T X aib (^^,e T \k\-,^,a F ,N S j .(17) 

Eq. (|16f) can then be written in terms of the dimensionless 
scaling variables s = hwf (fcgT), t = ^r|k| as 
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Therefore any solution of this equation must have the 
scaling form, s — s (t;—,a.F, N) , which is equivalent to 
Eq.©. For purely real X a j, this equation sets the bound 
for allowed frequencies < c|k|. 

To proceed, we need to resort to approximations 
for H a ,b- Consider first II a . Within the random 
phase approximation, we find that for s 2 > 

t 2 the polarization scaling function JSJ "P a (s,i) w 

~ln2t»-." /l_^J^ N l_ l VipE tailh! 

fTrTj) 



a F N 
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can be 



the above approximation, the Eq. 
solved by first ignoring the imaginary part of II a and 
then including it perturbatively. Then, to th or- 
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FIG. 1: (Black) The dispersion of the thermo-plasma po- 
lariton mode in terms of the dimensionless scaling vari- 
ables s = hu/k B T = Lj/(131GHz)(Kelvin/T) and t = 
hvF\k\/kBT = 6.1/im|k|(Kelvin/T). At small momenta, the 
dispersion asymptotes to c|k|- the light cone (red). Above the 
crossover scale t c w g(c—^/c? — v 2 F )/^c 2 — v 2 F , the dispersion 
is given by s = y/t(t+g)/y / t + 2g, where the effective electron- 
electron interaction strength is g = e 2 N In 2/ {^I-kJivf)- For 
orientation, at room temperature and for |k| _1 = 400nm, 
Lu p ~ 6THz. 



cxfN \n2/(2Tr). This is a cubic equation in s 2 and so 
a closed form solution exists. The expression is too long 
to present here but the solution is plotted in Fig^ In 
the low k regime the dispersion follows the light cone, 
i.e. in the scaling variables s = (c/vp)t. At a crossover 
scale t c w g(c - y 'c 2 - v 2 F )/ y 
changes to s = y/t(t + g)j\ft- 



the dispersion 
2g. The physical ori- 
gin behind these two regimes can be understood as fol- 
lows: if the interaction between the electrons were re- 
placed by an instantaneous Coulomb interaction, the 
plasma mode would continue dispersing as scouiomb = 
y/i(t + g)/y/t + 2g even for vanishing t. However, at the 
crossover scale t c , the light cone crosses scouiomb an d the 
two mix. This is the plasmon-polariton effect. 

As vf/c — > 0, the light cone regime regime shrinks as 
gyp /(2c 2 ) and 
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(19) 



The imaginary part of the frequency is the k— dependent 
inverse lifetime of the monochromatic plasmon-polariton 
mode and has its origin in the electron-hole scattering of 
the collective mode. The distance over which the mode 
propagates, the propagation length S£, has a scaling form 
as well, and it is given by the product of the lifetime and 
the phase velocity 81 = £t^s[s] / (t^sm[s]) . This quantity 
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FIG. 2: The thermo-plasma polariton (in-plane) propaga- 
tion length, S£, (dashed red) and the (out-of-plane) atten- 
uation length, £ a , (solid black) normalized to thermal length 
(It = hvF I '(feT)) vs. the mode frequency, oj p , normalized by 
ksT/h. At room temperature and for lo v = 20THz, 51 m 2/im. 



There certain similarities between the thermo-plasma 
polartitons discussed here and the surface plasma po- 
laritons in thin metals. In particular in both cases the 
evanescent electromagnetic fields prevent the power from 
radiating outwards. However, in the case of the surface 
plasma polaritons, the dispersion is essentially fixed by 
the electron density and does not depend significantly on 
temperature (or gate voltage). This is in sharp contrast 
to the strong temperature and gate voltage dependence 
of the properties of the thermo-plasma polaritons, the lo- 
cal manipulation of either of which may open new ways 
to generate and waveguide these modes. 
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is plotted in Fig. (J2J for N — 4 and vf — 8 x 10 5 m/s. 
Meanwhile, note that, within the same approximation for 
lib, there is no weakly damped solution to Ea. l(ro|) . 

Finally, given the scalar A and vector potential A4 
solution for the evanescent wave, the electric and mag- 
netic field distributions associated with the solution can 
be determined. While the magnetic field is always trans- 
verse, the electric field has a component perpendicular to 
the graphene sheet. Up to an overall scale, its in-plane 
Fourier components are 
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The fields are thus bound to the surface which prevents 
the power from radiating away from the graphene sheet. 

While most of the results presented here concern un- 
gated graphene, the effects of applied gate voltage can 
be easily taken into account within the scaling theory. 
Finite gate voltage V g introduces another energy scale 
which must be compared to fceT. If V g <C ksT then 
the results presented above hold. If V g » fc^T, then 
the thermal energy (length) scale fcgT, (#r), must be 
replaced by V g or by hvF/V g respectively. 

Experimentally, the temperature dependent attenua- 
tion of the EM fields perpendicular to the layer can be 
measured using near field optical microscopy, while the 
dispersion can be measured, for instance, by inelastic 
light scattering. 
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